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Abstract. Four XMM-Newton observations of the core of M31, spaced at 6 month intervals, show that 
the brightest point X-ray source, RX J0042.6+4115, has a 0.4-10 keV luminosity of ~5 10 38 erg s" 1 , 
and exhibits significant variability in intensity and X-ray spectrum over a time scale of ~100 s including 
hard flares; such behaviour is only observed in Z-sources and transient black hole binaries in our Galaxy. 
The lightcurves, X-ray spectra and hardness-intensity data from the four XMM-Newton observations all 
strongly suggest that it is a Z-source, bringing the total number of known Z-sources to nine. 
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1. Introduction 

The Andromeda galaxy (M31) is the nearest spira l 
galaxy, at a distance of 760 kpc l|van den Bergrl20 00) . 
The X-ray emission from M31 is dominated by point 
sources. While the study of variability in X-ray sources 
in external galaxies has been limited by the sensitiv- 
ity of observatories prior to XMM-Newton, its three 
large X-ray telescopes can make variability within a 
single observation visible down to fluxes of ~10~ 14 er 



cm 



in the 0.3 



)le down to fluxes oi ~iu erg 
10 keV band (e.g. Ijansen et al.l 



l200llL equivalent to luminosities of ~10 erg 
M31. The XMM-Newton observations of M31 have 
already revealed a persistent, stellar-mass black hole 
binary l|Barnard et alJ l2003al Paper 1), periodic in- 
tensity dips in the M31 globular cluster X-ray source 
Bol58 (jTrudolvubov et al.ll2002l) and pulsa tions in a 
white-dwarf transient l|Osborne et al.ll200l[) . 

RX J0042. 6+4115, the subject of this pa- 
per, was discovered in th e Einstein survey of 
Trinchieri fc Fabbian 3 <ll99lh . and was named by 



Supper et al.1 1) 1 997(1 . who associated it with a globu- 
lar cluster. However, its location w as more accurately 
given in a recent Chandra survey l|Kong et alJ 12002) 
as 00 h 42 m 38*581 +41° 16' 03'.'80 with a positional 
error of O'.'Ol, and was not associated with any fore- 
ground object, glo bular cluster or background AGN. 
iKong et alJ {2002) found variability in the lightcurves 
and spectra of RX J0042. 6+4115 and report a lumi- 
nosity in the 0.3-7 keV band of 1.5 1 38 erg s" 1 , as- 
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suming a standard absorbed power law spectral model 
with an index of 1.7 and absorption equivalent to 10 21 
hydrogen atom cm -2 . Our results from four XMM- 
Newton observations of RX J0042. 6+4115 reveal 0.3- 
7 keV luminositie s of ~3-4 times that reported by 
IKong et alJ l|2002h while spectral fits required a black- 
body component in addition to a power law compo- 
nent. Furthermore, significant variability is observed 
over time-scales of ^100 s. Of the Galactic X-ray bi- 
nary population, such behaviour is only observed in Z- 
sources and X-ray transients. As this object is clearly 
not a transient, we argue that RX J0042. 6+4115 is the 
9th Z-source. 

The paper is constructed as follows: the prop- 
erties of Z-sources are briefly reviewed in the next 
section; the XMM-Newton observations and analy- 
sis techniques are then discussed; lightcurves, colour- 
colour diagrams and energy spectra are presented 
in the Results section and the case for classifying 
RX J0042. 6+4116 as an analogue to the Galactic Z- 
sources presented in the Discussion section. 



2. Properties of Z-sources 

The Z-source subclass of low mass X-ray binaries 
(LMXB) is ch aracterised by high X-ray lum inosities 
(>10 38 erg s~ 1 . lHasinger fc van der Klisll989h and in- 
tensity v ariations of up to 100% over timescales of 
minutes ijWhite et alJ Il97fl and references within); 
they are named after the patterns traced on their 
colour-colour diagrams. If X-ray lightcurves of a 
source are obtained in several energy bands, then 
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Table 1. Journal of XMM-Newton observations of the 
M31 core 



Fig. 1. A detail of the 0.3-10 keV PN image from 
the June 2001 XMM-Newton observation of the core 
of M31; north is up, east is left. The image is log scaled 
and 5' across. RX J0042. 6+4115 is circled in black, the 
circle defining the source extraction region. 



one can define hard and soft colours to see spec- 
tral variations over time; the commonly used colour- 
colour diagram (CD) is a plot of hard colour vs. 
soft colour. Z-sources exhibit three-branched, approx- 
imately Z-shaped CDs, with temporal properties that 
are stro ngly correlated to the position o f the source on 
the CD l|Hasinger fe van der Klisll 989). The branches 
are historically named the horizontal branch (HB), 
normal branch (NB) and flaring branch (FB). During 
spectral evolution, the Z-source moves smoothly along 
the branches without jumping from one branch to 
another; the mass accretion rate is thought to in- 
crease as the source travels in the s ense HB to NB 
to FB l|Hasinger fc van der Klisl ll9891. The Z-sources 
exhibit significant branch movement over time-scales 
of ^100 s or les s and trace out th eir complete pat- 
tern in -1 day l)Muno et alJl2002t . All Galactic Z- 
sour ces exhibit kiloh ertz quasi-periodic oscillations 
(see iLiu et ahl l200li and references within). This 
tells us that they all contain neutron star primaries, 
since the power density spectra of LMXB with black 
hole primarie s are not thought to contain any power 
at >500 Hz ijSunvaev fc Revnivtsevll2000l) . These Z- 
sources may be divided further into two subsets: 
CygnusX-2, GX5-1 and GX 340+0 have strong, hor- 
izontal HBs and weak flaring branches while ScoX-1, 
GX17+2 an d GX 349+2 have strong FBs a nd weak, 
vertical HBs l|Hasinger fc van der KliJl 989b The un- 
derlying differences that account for this division of 
Z-sources is unknown. Two additional Z-source candi- 
dates have recen tly been identified: Circinus X-l (e.g. 
Ilaria et aL |200l| ) and the first extr a-galactic Z-source 
LMCX-2 Jsmale fc Kuulkenfcood) . It is highly likely 
that the set of 6-7 known Z-sources in our Galaxy 
is complete, due to their high X-ray luminosities, yet 
they account fo r only ~ 5% of the Galactic LMXB 
I Liu et alJEoOljl . 



Observation 


Date 


Exp 


Filter 


1 


25/07/00 (revOlOO) 


34 ks 


Medium 


2 


27/12/00 (rev0193) 


13 ks 


Medium 


3 


29/06/01 (rev0285) 


56 ks 


Medium 


4 


06/01/02 (rev0381) 


61 ks 


Thin 



The spectral evoluti on of ScoX-1 a l ong its flaring 
branch was studied bv iBarnard et al.l ll2p03bh. using 
data from the Rossi-XTE observatory l|Bradt et al.l 
Il993h : luminosities of ~5-12 1 38 erg s 1 were ob- 
served in the 1-30 keV energy band. This exceeds the 
Eddington luminosity for an accreting, 1.4 M Q neu- 
tron star, the luminosity where the radiation pres- 
sure on accreting mater ial balances the gravitation 
forces (~2 10 38 erg s- 1 ), iFrank et alJ 2002L However, 
the majority of the X-ray emission comes from an 
extended accretion disc corona and only 10-30% of 
the 1-30 keV emission comes from the neutron star 
surface, in the form of blackbody emission; the lumi- 
nosity of the neutron star is generally sub-Eddington 
l|Barnard et al.ll2003bj) . We can therefore expect the 
X-ray spectrum of an analogue of ScoX-1 to contain 
a strong blackbody component as well as the stand ard 
power law component fitted bv lKong et all l|2002h . 



3. Observations 

Four XMM-Newton observations were made of the 
core of M31; data from the three EPIC detectors 
in each observation were analysed for variation be- 
tween observations and variability within observa- 
tions. The three im aging CCD in struments, two MPS 
l|Turner et alJl200lJ) and one PN Intruder et al.ll200ll) 
share a common, circular field of view with a radius 
of 30' l)jansen et alJl200llh A journal of the XMM- 
Newton observations of the core of M31 are presented 
in Tabled Lightcurves of RXJ0042. 6+4115 were ob- 
tained in the 0.3-2.5, 2.5-10, 4.0-10 and 0.3-10 keV 
energy bands from the PN and MOS detectors at their 
highest time resolutions. The analysis closely follows 
that performed in Paper 1, with two minor differences. 
Firstly, the background contribution to the lightcurves 
was very small, and all good data were used. Secondly, 
spectra from the PN only were analyzed; the data were 
accumulated into 250 eV bins and bins below 0.4 and 
above 10 keV were ignored. 

4. Results 

Observations 1, 3 and 4 exhibit significant variabil- 
ity in all energy bands, while the lightcurves from 
observation 2 exhibit no significant variability; each 
lightcurve was fitted with a line of constant inten- 
sity and the resulting best fit x 2 /dof and probabil- 
ity of randomly getting a larger x 2 /dof (P) are pre- 
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Table 2. Variability of EPIC lightcurves of RX J0042.6+4115 in different energy bands with 400 s binning: 
X 2 /dof and P (the probability that a random fit to the data would be a worse fit) are shown for the best fit 
lines of constant intensity to lightcurves in three energy bands for each observation. 



Observation 


Xo.3-2.5kcv/dof 


P().3-2.5kcV 


X2.5-10kov/ dof 


P2.5-10kcV 


Xo.3-iokcv/dof 


P().3-10keV 


1 


100/77 


0.04 


143/77 


7 10" 6 


162/77 


6 10" 6 


2 


33/24 


0.10 


23/24 


0.52 


28/24 


0.26 


3 


513/131 


<HT 6 


412/131 


■C 10~ 6 


259/131 


<10~ 6 


4 


730/152 


<10" 6 


497/152 


<10~ 6 


441/152 


<10" 6 




10* 2x10* 3x10* 4x10* 5x10* 6x10* 

Time (s) 

Fig. 2. The EPIC (PN + MOSl +MOS2) lightcurves 
of RXJ0042.6+4115 from observation 4 in the 0.3- 
2.5 keV (top panel) and 2.5-10 keV (bottom panel) 
energy bands; the upper curve in each panel is the 
background-subtracted source curve, while the lower 
curve is the background curve. The lightcurves were 
accumulated with time bins of 400 s. Clearly, the am- 
plitude of flaring is higher in the 2.5-10 keV energy 
band, and it is also clear that the background makes 
no contribution to the flares. 



sented in Table Strong intensity enhancements oc- 
cur at the start of observations 3 and 4; in fact, the 
lightcurves of RX J0042. 6+4115 closely resemble the 
1998, January RXTE lightcurves of ScoX-1, which 
captured a transition between the FB and the NB 
l)Barnard et alJl2003bj) . The flares in ScoX-1 last for a 
few thousand seconds and are hard, meaning that the 
fractional amplitud e of the flares increases with energy 
llWhite et al.lll976h . Combined EPIC (PN + MOSl + 
MOS2) lightcurves of RX J0042. 6+4115 from observa- 
tion 4 in the 0.3-2.5 and 2.5-10 keV energy bands are 
given in Fig. with 400 s binning; clearly the relative 
amplitude of flaring is higher in the 2.5-10 keV energy 
band (~100%) than in the 0.3-2.5 keV band (~ 33%). 

Hardness-intensity data was obtained for the 4 ob- 
servations of RX J0042.6+4115; the hardness was de- 
fined as the ratio of the 4.0-10 keV intensity to the 
0.3-2.5 keV intensity, and data were integrated over 
2000 s. The high energy band 4.0-10 keV was cho- 
sen to better separate the three branches; Z-sources 
are more varia ble at higher energ ies than at low en- 
ergies (see e.g. IWhite et alJfl97fil who modelled the 
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Fig. 3. Hardness ratio vs. intensity for 4 XMM- 
Newton observations of RX J0042. 6+4115; the hard- 
ness ratio is defined as l4.o-iokcv/Io.3-2.5kcV and the 
data is binned to 2000 s. The dashed line shows a pos- 
sible Z-track for RX J0042. 6+4115. 



X-ray spectrum as bremsstrahhmg emission and in- 
terpreted the increase in hardness with intensity as 
an increase in temperature). Fig. shows the hard- 
ness vs intensity diagram (HID) for all observations of 
RX J0042. 6+4115; stars, circles, triangles and squares 
represent data from observations 1, 2, 3 and 4 respec- 
tively. The resulting HID appears to show fragments 
of the Z-p attern characteristic of the Cvg nusX-2 like 
Z-sources l Hasinger fc van der Klidll98fl|) : it appears 
that observation 1 is on the NB, observation 2 is near 
the HB/NB apex, observation 3 is on the FB, and 
observation 4 is on the HB. 

In order to illustrate the effects of long integration 
times on HID for Z-sources, we present in Fig. ^hard- 
ness vs. intensity data from 400 ks of RXTE data on 
ScoX-1 using 96 s binning and 2000 s binning. We see 
that using 2000 s bins preserves the shape of the HID 
that would be seen at 96 s resolution, although the 
curve is more sparsely sampled; hence we can be con- 
fident that our HID of RX J0042. 6+4115 is consistent 
with Z-source HID. 

The intervals of high intensity variability were 
omitted from spectral analysis, because the hardness 
ratios of these data (and therefore spectra) signifi- 
cantly varied; as a result, the spectrum from obser- 
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2.5x10* 3xl0 4 
Intensity (count s -1 ) 

Fig. 4. Hardness vs. intensity data from RXTE obser- 
vations of the Galactic Z-source ScoX-1 with 96 s bin- 
ning (grey crosses) and 2000 s binning (circles). The 
HID at 2000 s binning clearly reflects the branch move- 
ment, even though the HID at 96s demonstrates that 
spectral variability occurs over much shorter time- 
scales. 



vation 1 had an exposure of only 10 ks, while the 
spectra of observations 3 and 4 both had exposures 
of —50 ks. The PN spectra from the four observations 
of RX J0042. 6+4115 were first fitted by simple power 
law models (A1-A4). None of these fits were statisti- 
cally acceptable. Fits to the spectrum of observation 
4 using models consisting of a single bremsstrahlung 
component (x 2 /dof=103/29) or a single COMPTT com- 
ponent (x 2 /dof = 247/27) were also rejected. Hence a 
two component emission model was applied, consist- 
ing of a blackbody and a power law. These two compo- 
nent models (B1-B4) all provided good fits to their re- 
spective spectra. All models incorporated freely-fitted, 
line-of-sight absorption. The unfolded, 0.4-10 keV PN 
spectrum of RX J0042. 6+4115 from observation 1 is 
presented in Fig. [SJ the data are shown (filled cir- 
cles) along with the power law (1) and blackbody (2) 
components, and the sum for the best fit model. The 
best fit parameters, x 2 /dof and fluxes of each model 
are presented in Table 90% confidence errors were 
obtained for the two-component models, but not the 
single component models since the fits were unaccept- 
able. The spectra were then tested for the presence of 
line emission. There is evidence for a 6.4 keV line, cor- 
responding to fluorescent iron Ka emission in obser- 
vation 3 only: F-testing showed that the improvement 
in x 2 /dof due to fitting the line had a 3% probability 
of being random. 

The spectral shape parameters from the two com- 
ponent models are fairly consistent from observa- 
tion to observation, with line-of-sight absorption of 
-0.14 1 22 H atom cm 2 , a blackbody tempera- 
ture equivalent to kT —1.1 keV, and a photon in- 
dex of —2.1. However, the normalisations of the two 
emission components change significantly. This is en- 



0.5 1 2 5 10 

Energy (keV) 

Fig. 5. The PN spectrum of RX J0042. 3+4115 from 
observation 1; the data (filled circles) are fitted with 
the sum of two components, a power law component 
(1) and a blackbody component (2). 

tirely consistent with t he spectral evolution of ScoX-1 
^Barnard et all2003b|) . The total 0.4-10 keV luminos- 
ity ranges over 4.8-5.6 10 38 erg s _1 , or —2.5 times the 
Eddington limit for a 1.4 M© neutron star, while the 
blackbody contribution is 0.7-1.6 1 38 erg s , which 
is sub-Eddington. The total 0.4-10 keV luminosity de- 
creases by —10% over the 4 observations, while the 
blackbody luminosity increases by 100%. The pecu- 
liar softness of observation 3 data appears to be due 
to the steep photon index, bearing out the associa- 
tion of observation 3 data with the flaring branch - 
the phot on index steepens alon g the flaring branch in 
ScoX-1 (|Barnard et alJl2003b|) . Meanwhile the hard- 
ness of observation 4 data may be explained by the 
increased blackbody contribution to the X-ray flux, 
since its fractional contribution to the 4.0-10 keV flux 
is greater than for the 0.3-2.5 keV flux, as demon- 
strated in Fig. 

5. Discussion 

In order to place RX J0042. 6+4115 in M31, we must 
show that it is not a foreground object or a background 
AGN. If RXJ0042. 6+4115 was local, it would have to 
be within a few kpc, and so the 0.4-10 keV luminos- 
ity would be reduced by a factor of — 10 6 , to —5 10 32 
erg s -1 . Thus RX J0042. 6+4115 would be too faint 
to be a persistently bright LMXB (|Lewin et alJll995h 
and would be either a cataclysmic variable or soft X- 
ray transient (SXT) in quiescence. However, the ob- 
served X-ray spectra for RX J0042. 6+4115 are com- 
plet ely unlike any observ ed in cataclysmic variables 
(see iKuulkers et all2003l and references within) . Also 
the blackbody temperatures in our spectra are many 
times high er than expected for a quies cent SXT (c.f. 
~ 0.3 keV, iMenou fc McClintockll200l|) . Finally, they 
are inconsistent with observed spectra from AGN, 
which in the 0.5-10 keV energy range are described 
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Table 3. Results from fitting the 0.3-10 keV PN X-ray spectra of RX J0042.6+4115 with spectral models A[n] 
and B[n]; A1-A4 are single component power law models fitted to observations 1-4 respectively, and B1-B4 are 
equivalent two-component models, consisting of a blackbody and a power law. Best fit parameters are given, 
along with the x 2 /dof and unabsorbed 0.4-10 keV fluxes of the total emission and the blackbody component, 
where applicable. Numbers in parentheses indicate 90% confidence errors on the last digit. 



Model 


at / 1 ri22 

A H /10 


kT 


ubb /10 


r 


/in— 3 

npo/10 


2 /J r 

X /dot 


rntot /in — 12 

b /10 


b /10 




atom cm -2 


keV 


erg s~ 
at 10 kpc 




photon cm -2 
at 1 keV 




erg cm -2 s _1 


erg cm -2 s _1 


Al 


0.14 






1.9 


1.4 


49/29 


8.1 




A2 


0.15 






1.8 


1.4 


49/24 


8.0 




A3 


0.14 






1.8 


1.2 


351/30 


7.3 




A4 


0.14 






1.8 


1.16 


258/29 


7.2 




Bl 


0.14(3) 


1.08(10) 


1.28(18) 


2.02(3) 


1.36(3) 


34/27 


8(3) 


1.1(4) 


B2 


0.14(2) 


1.05(8) 


1.87(15) 


1.95(4) 


1.22(3) 


23/22 


8(2) 


1.5(5) 


B3 


0.16(2) 


1.11(3) 


2.79(8) 


2.24(2) 


1.12(1) 


37/28 


7.3(8) 


2.3(3) 


B4 


0.13(2) 


1.08(3) 


2.65(10) 


2.10(3) 


1.00(2) 


33/27 


7.0(9) 


2.2(3) 



by a simple power law l|Pagelll99 8) whereas the spec- 
tra of RX J0042.6+4115 require a second component. 
It would be hard to understand this object if it were in 
our Galaxy or in a distant background galaxy, whereas 
it is easily comprehensible as a member of M31. 

RX J0042. 6+4115 is therefore a high luminosity X- 
ray binary, with a total 0.4-10 keV luminosity of 5- 
6 10 38 erg s -1 , including 0.8-1.7 1 38 erg s _1 from 
a ~1 keV blackbody component. The total luminos- 
ity exceeds the Eddington limit for a 1.4 M Q neutron 
star, however the blackbody component does not. Two 
Galactic Z-sources are known to regularly e xceed the 
Eddi ngton limit: ScoX-1 and GX5— 1 (e.g. iLiu et al.l 
2001), w hile of the black h ole XB, only transients in 



outburst I Chen et alJll997ft or extreme odd balls such 
as GRS 1915+105 (e.g. iBelloni et al.ll2000h reach the 
luminosities we have observed in RXJ0042. 6+4115. 

Outbursts in X-ray transients generally rise in in- 
tensity quickly (a few days) and decay exponentially; 
the decline time-scales cluster around 30 days, al- 
thoug h time-scales of 400 days are known llChen et al.l 
119971) . The first detection of RX J0042. 6+4115 was 
over 20 years ago by the Einstein satellite, when 
the 0.2-4.0 keV luminosity was 2.4 1 38 erg s^ 1 



[ Trinchieri fe Fabbianol Il99l|h and it has been per- 
sistently bright in ev ery observation since (e.g. 
ISunner et al.lll997lEoOlh . and so is highly unlikely to 
be a transient. 

Since the three branches of the Z were observed 
over 18 months, it is possible that RXJ 0042. 6+4115 
is a ve ry luminous atoll source. However. iMuno et al.l 
l|2002() report that three-branched colour diagrams 
in atoll sources are only observed in systems where 
the intensity varies by a factor of >80; atoll sources 
that vary in X-ray intensity by no more than a fac- 
tor o f 10 only ever d isplay a portion of the pat- 
tern l|Muno et alJl2002|) . The 0.3-10 keV intensity of 
RXJ0042.6+4115 varies by a factor of -2 over the 
four XMM-Newton observations, and yet three pat- 
terns of hardness-intensity behaviour are observed. In 



this regard, RX J0042. 6+4115 is unlike any Galactic 
atoll source. 

The lightcurves, hardness- intensity diagram, and 
the luminosities and composition of the X-ray spec- 
tra all strongly suggest that RX J0042. 6+4115 is a Z- 
source; its luminosity and Z-pattern suggest that it is a 
super-Eddington, CygX-2-like Z-source analogous to 
GX5-1. If RXJ0042.6+4115 is a Z-source, it would 
be only the 9th to be found. 

Acknowledgements. We thank the referee for providing 
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supported by PPARC. Data from the HEASARC public 
data archive for RXTE were used in constructing Fig. fj] 
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